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I .  SUMMARY 

This  report  describes  the  operation,  maintenance  and  research  activities  .it 
the  Norwegian  Seismic  Array  (NORSAR)  for  the  period  1  April  to  30  September 
1983. 

The  uptime  of  the  NORSAR  online  detection  processor  system  has  averaged  9S.H 
as  compared  to  95.7  for  the  previous  period.  Most  of  the  downtime  was  caused 
consisted  of  short  breaks  caused  by  resync  of  lines  and  correction  of  TO  Li- 
caused  M0DC0MP  restart.  A  total  of  2295  events  were  reported  in  this  period, 
giving  a  dalLy  average  of  L2.5  events.  The  number  of  reported  events  per 
month  varies  from  264  in  April  to  463  in  May.  There  have  been  some 
difficulties  with  the  communications  lines;  03C  was  affected  several 
times  by  power  line  damage,  04C  was  down  during  part  of  the  period  due 
to  a  damaged  cahle,  and  02C  was  down  the  last  part  of  the  period  due  to 
problems  on  the  communications  line. 

Improvements  and  modifications  in  the  NORSAR  on-line  system  are  briefly 
described  in  Section  HI.  Quality  and  uptime  for  the  plotting  system  has 
been  Improved  by  use  of  a  Versatec  raster  plotter  instead  of  OaLr.omp  plot!  " 
A  new  disk  drive  lias  also  been  installed  to  increase  the  on-line  data 
capacity  from  30  to  45  hour'-.  The  ARPANET  connection  is  again  available 
from  NORSAR. 

Section  IV  describes  field  insr rumentat ion  and  maintenance  activities  a‘ 
the  NORSAR  Maintenance  Center  and  Includes  an  overview  of  the  status  of 
the  NORESS  and  02B  telemetry  stations 

The  research  activity  is  briefly  described  in  Section  VI.  Subsection  i 
presents  an  example  of  a  seismic  absorption  band  at  high  f reqene ios .  Sub¬ 
section  2  discusses  spectra'  bandwidth,  pulse  width  and  moments  in  source 
analysis.  Subsection  3  gives  a  description  of  experiments  with  temporary 
field  installations  in  the  new  regional  array  during  the  summer  of  1983 
and  some  preliminary  results  from  analysis  of  data  collected.  Subsection 
4  discusses  further  developments  In  the  Regional  On-line  Array  ro-ossf  tig 
Package  (RONAPP).  Experiments  Involving  weighted  bean  form  lag  in  t  real 
time  environment  are  presented  in  subsection  5. 


II.  OPERATION  OF  AI.L  SYSTEMS 
a .  1  Detection  Processor  (DP).  Operation 

There  Have  been  138  breaks  in  the  otherwise  continuous  operation  of  the 
NOKSAP  online  system  within  the  current  6-month  reporting  interval.  60% 
of  the  stops  were  short  breaks  that  lasted  for  no  more  than  two  minutes. 

These  breaks  were  not  caused  by  the  complete  DP  system  going  down,  rather 
by  resync  of  lines  and  correction  of  TOD-eaused  M0DC0MP  restart.  The 
uptime  percentage  for  the  period  is  98.8  per  cent  as  compared  to  93.7 

for  the  previous  period. 

Fig  11.1.1  and  the  accompanying  Table  II. 1. 1  both  show  the  daily  DP  downtime 
for  the  days  between  1  April  and  30  September  1983.  The  monthly  recording 
tLnes  and  percentages  are  given  in  Table  II. 1.2. 

The  breaks  can  be  grouped  as  follows: 

a)  Stops  related  to  program  work  and  errors 

2 

b)  Hardware  maintenance 

c)  Power  breaks 

33 

d)  TOD  correction 

e)  Array  monitoring 

*  ,,  14 

f)  M0DC0MP  failure 

g)  Disk  failure  1 

The  total  downtime  for  the  period  was  54  hours  and  36  minutes.  The  mean- 
tlme-between-fai lures  (MTBF)  was  1.3  days  as  compared  with  2.9  days  for  the 

previous  period. 
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Month 

DP 

DP 

No.  of 

No .  of 

DP 

Uptime 

Upt ime 

DP  Breaks 

Days  with 

MTBF* 

(hrs) 

(%) 

Breaks 

(days) 

Apr 

697.15 

96.8 

32 

9 

0.9 

May 

737.00 

99.1 

9 

6 

3.1 

Jun 

702.03 

97.5 

11 

9 

2.4 

.Jul 

743.42 

99.9 

14 

10 

2.1 

Aug 

739.25 

99.4 

46 

15 

0.7 

Sep 

713.46 

99.8 

26 

12 

1.1 

4338.11 

98.8 

138 

61 

1.3 

♦Mean- 

time-betueen- 

failures  = 

(Total  uptime/No.  of  up 

intervals) 

TABLE  11.1.2 

Online  System  Performance 
1  April  -  30  September  1983 


II. 2  Event  Processor  Operation 

In  Table  II. 2.1  some  monthly  statistics  of  the  Event  Processor  operation  are 
given: 


Teleseismic 

Core  Phases 

Sum 

Daily 

Apr  83 

264 

60 

324 

10.8 

May  83 

463 

56 

519 

16.7 

Jun  83 

317 

67 

384 

12.8 

Jul  83 

313 

82 

395 

12.7 

Aug  83 

320 

56 

376 

12.1 

Sep  83 

265 

32 

297 

9.9 

1942 

353 

2295 

12.5 

TABLE  II. 2.1 


B.  Kr.  Hokland 
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II. 3  Array  Communication 

Table  II. 3.1  reflects  the  performance  of  the  communications  system  through¬ 
out  the  reporting  period. 

High  figures  as  indicated  in  the  table  are  related  to  incidents  such 
as  power  line  damage  (03C),  NORESS  usage  of  the  original  06C  line,  and 
damaged  cable  (04C).  Otherwise  communication  systems  have  been  most 
reliable. 


Summary 

April :  Apart  from  an  incident  in  connection  with  06C  communication 

system  (week  16),  reliable  performance. 


May:  NORESS  data  transferred  via  the  06C  original  communication 

path;  otherwise  reliable  performance. 

June :  Unwanted  noise  observed  on  the  04C  communication  line  (27  June). 

Test  carried  out  (29  June)  did  not  revea.  irregularities,  but 
situation  changed  (30  June).  Transmitted  commands  had  no  in¬ 
fluence  on  the  CTV  modem.  NTA/Lillestrom  notified.  06C  comnunica- 
tlon  path  still  in  NORESS  usage. 

July :  04C  back  in  operation  after  NTA/Haraar  had  replaced  part  of  the 

communication  cable  between  the  CTV  and  the  local  central. 


August :  A  line  problem  caused  02C  outage  between  12-lb  August.  03C 

dropped  out  (29  August),  and  remained  so  for  the  rest  of  the 
period.  Remaining  system  reliable  performance. 

September :  03C  back  in  operation  15  September  after  a  power  outage.  The  same 
subarray  lost  its  power  again  (26  September,  also  this  time  caused 
by  falling  trees  over  the  power  line).  02B  was  out  of  operation 
27-28  September  due  to  NTA/ Lillestrtfra  work. 
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Table  II. 3. 2  indicates  distribution  of  outages  with  respect  to  the  individual 
subarrays. 

Miscellaneous 

The  new  line  to  be  used  in  connection  with  the  NORESS  array  was  not  opera¬ 
tional  until  mid-July,  although  NTA/Hamar  and  Lillestrftn  were  engaged  in 
May  and  June.  Main  reason  for  the  delay  was  unexpected  attenuations  in 
parts  of  the  line  between  Hamar  and  the  array. 

As  far  as  the  temporary  NORESS  array  and  02B  (telemetry)  installation  is 
concerned,  we  have  observed  a  few  minor  problems,  such  as:  spikes  and 
reduced  performance  on  NORESS  channels  02,  03,  12  and  13  (week  19) 

Reduced  performance  on  ch.  12  (NORESS)  (week  39). 

NORESS  channels  17-22  out  of  operation  for  2  days  (week  39).  Heavy 
attenuation  on  ch.  27,  02B  (telemetry).  Sept  1  (week  39). 


O.A.  Hansen 
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Sub- 

Array 

APR  (4) 
(4. 4-1. 5) 

MAY  (5) 
(1.5-5. 6) 

JUN  (4) 

(6. 6-3. 7) 

JUL  (4) 
(4-31.7) 

AUG  (5) 
(1.8-4. 9) 

SEP  (4) 
(4.9-2.10) 

AVERAGE 
i  YEAR 

01A 

0.005 

0.021 

0.028 

0.001 

0.005 

0.005 

0.011 

01B 

0.004 

0.023 

0.027 

0.002 

0.005 

0.005 

0.011 

02B 

0.005 

0.023 

0.027 

0.002 

0.004 

3.58 

0.61 

02C 

0.005 

- 

0.093 

0.002 

0.71 

0.004 

0.14 

03C 

0.005 

0.055 

0.128 

*3.52 

*20.0 

*52.7 

*12.73 

04C 

0.005 

0.022 

*7.714 

*33.08 

0.14 

0.01 

*6.87 

06C 

0.23 

*2.07 

*100.0 

*81.09 

*100.0 

*100.0 

*63.89 

AVER 

0.04 

0.3 

15.4 

16.8 

17.3 

22.3 

12.0 

04C 

03C 

03C 

03C 

LESS 

06C 

06C 

06C 

06C 

06C 

- 

- 

1.3 

0.7 

0.2 

0.7 

1.5 

*  See 

item  11.3 

(array  communications) 

regarding 

figures  with  asterisks 

TABLE  XI. 4.1 

Communications  performance.  Figures  in  per  cent,  based  on 
total  transmitted  frames/week  (1.2096  x  10^). 

(4  Apr  -  30  Sep  83) 


5 

3 
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Week/  Subarray/per  cent  outage 


Year 

01A 

01B 

02B 

02C 

03C 

04  C 

06  C 

14/83 

0.020 

0.020 

0.020 

0.020 

0.020 

0.020 

0.020 

15 

- 

- 

- 

- 

- 

- 

- 

16 

- 

- 

- 

- 

- 

- 

4.3 

17 

0.020 

0.001 

0.002 

- 

0.002 

0.001 

0.9 

18 

- 

- 

- 

- 

- 

- 

6.4 

19 

- 

- 

- 

- 

- 

- 

3.55 

20 

0.070 

0.060 

0.090 

- 

0.015 

0.060 

0.4 

21 

0.035 

0.020 

0.025 

- 

0.120 

0.050 

0.04 

22 

- 

- 

- 

- 

0.04 

0.03 

100.0 

23 

0.08 

0.11 

0.074 

0.370 

0.085 

0.08 

100.0 

24 

0.03 

0.002 

0.04 

- 

0.04 

0.03 

100.0 

25 

- 

- 

- 

- 

0.38 

0.13 

100.0 

26 

- 

- 

- 

- 

0.01 

30.6 

100.0 

27 

0.002 

0.0023 

0.0024 

0.0015 

0.002 

71.4 

100.0 

28 

- 

- 

- 

- 

0.08 

60.8 

76.1 

29 

0.0001 

0.0013 

0.001 

0.0013 

0.0013 

0.005 

48.3 

30 

0.002 

0.003 

0.003 

0.0034 

0.004 

0.10 

100.0 

31 

0.0002 

- 

- 

0.0004 

- 

0.014 

100.0 

32 

- 

- 

- 

- 

- 

0.2 

100.0 

33 

0.001 

0.0015 

0.001 

3.6 

0.001 

1.8 

100.0 

34 

0.001 

0.001 

0.001 

0.0014 

0.003 

0.03 

100.0 

TABLE  II. 4. 2 
Subarray/percent  outage 
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III.  IMPROVEMENTS  AND  MODIFICATIONS 

III.l  NORSAR  on-line  system  using  IBM  4531/4341  and  MODCOMP  Classic 
We  refer  to  the  detection  processor  operation  statistics  for  detailed 
performance  of  the  new  system. 

During  this  reporting  period  there  have  been  unit  checks  between  MODCOMP 
and  the  4331  computer.  The  follow-up  result  of  this  is  periodic  loss 
of  a  0.5  sec  data  block.  At  the  moment  we  cannot  say  where  the  problem 
is  located.  However,  we  expect  a  solution  of  the  communication  problem 
when  we  interface  the  MODCOMP  directly  to  the  4331  channel.  At  the 
moment  we  have  connected  the  MODCOMP  computer  directly  to  the  channel 
and  the  operating  system  on  IBM  has  accepted  the  device  which  simulates 
a  tape  drive.  We  have  also  done  some  read/write  tests. 

The  quality  and  uptime  for  the  plotting  system  has  been  improved.  We  are 
now  using  a  10"  Versatec  raster  plotter  instaed  of  the  Calcomp  plotters. 
The  Versatec  is  connected  through  the  Series/1  front-end  processor. 

The  4341  system  does  the  vector-to-raster  converting  and  transmit  raster 
lines  to  the  Series/1  processor  which  does  the  plotting  on  the  Versatec. 

The  ARPANET  connection  is  again  available  from  NORSAR.  At  the  moment 
we  have  a  direct  line  to  a  VAX-computer ,  which  Is  connected  to  the 
ARPANET  as  a  host  machine.  The  user-id  for  NORSAR  is  FRODE  and  the  machine 
id  is  NTA-VAX.  The  next  step  in  this  ARPANET  connection  is  a  host  con¬ 
nection  at  NORSAR.  We  will  probably  use  an  IBM  personal  computer  for 
this.  When  we  have  the  host  connection,  it  will  be  possible  to  transfer 
data  files  directly  from  the  main  system  through  the  ARPANET. 

We  have  also  Installed  a  new  disk  drive  to  keep  more  data  on-line. 

At  the  moment  we  have  a  capacity  of  30  hours  with  on-line  data.  With 
the  new  disk  drive  we  may  store  45  hours  with  on-line  data. 


R.  Paulsen 
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IV.  FIELD  INSTRUMENTATION  AND  MAINTENANCE  ACTIVITIES 
Improvements  and  modifications 

Reference  Is  made  to  Table  IV. 1  and  IV. 2  indicating  the  status  of  the  SP 
instruments  (original  array),  the  modified  NORESS  array,  and  the  expanded 
02B  subarray  (telemetry  stations). 

Since  the  last  report  the  NORESS  array  has  been  subject  to  a  number  of 
changes. 

5  June:  Five  sites  (see  Table  VI. 2)  were  equipped  with  3-component  SP 
seismometers  S-13,  as  a  temporary  experiment  to  collect  SP  data  recorded 
at  NORSAR.  The  array  was  operational  10  June,  and  the  operation  was 
terminated  3  July.  Otherwise  refer  to  Fig.  IV. 1  and  Table  IV. 3,  which 
indicate  station  numbering  and  geographical  coordinates,  respectively. 

23  July:  The  new  temporary  NORESS  array  was  fully  operational,  equipped 
with  vertical  seismometers  S-13  at  12.5  Hz,  occupying  NORSAR  analog  channels 
1-21,  wind  speed  data  on  ch.  22  (see  Table  IV. 2).  Reference  is  also  made 
to  Fig.  IV. 2,  and  Table  IV. 4,  which  indicate  station  numbering  and 
geographical  coordinates,  respectively. 


02B  (telemetry)  st.  1-6  equipped  with  seismometers  S-500,  now  occupy 
ch.  23-28  (see  Table  IV. 2  and  in  addition  Fig.  IV. 3  and  Table  IV. 5). 

Planning  with  regard  to  the  permanent  NORESS  is  going  on.  Meetings  have 
been  held  between  representatives  from  Sandia,  a  consultant  company 
and  NORSAR. 

Activities  in  the  field  and  at  the  Maintenance  Center 
This  section  outlines  in  brief  the  activities  in  the  field  and  at  the 
NORSAR  Maintenance  Center.  Table  IV. 6  indicates  other  activities  in  the 
field  and  at  the  NMC. 


Corrective/preventive  maintenance  in  the  original  array  has  been  limited 
to  replacement  of  33  RA-5  amplifiers  in  the  well  head  vaults,  and  replace¬ 
ment  of  a  remote  centering  device  (mass  position)  on  an  LP  instrument. 


Array  status 

As  of  30  September  1983  the  following  channels  deviated  from  tolerances: 
01B05;  02B05;  02C06;  03001,08;  06C  all  channels  (not  used). 


01A  01  8  Hz  filter 

02  ,  60  m  hole 

03  Wind  speed  measurements 

04  Attenuated  40  dB 

05  Wind  direction  measurements. 


O.A.  Hansen 


l 


Subarray 

(Normally) 

Instr.  no. 
within  SA 

Ch .  no .  on 
NORSAR  data 
tape 

Status  /  Time  of  changes 

OiA 

1 

1 

8  Hz  filter  (10/29/82) 

(1) 

2 

2 

1) 

(10/08/82) 

3 

3 

2) 

(10/28/82) 

4 

4 

5 

5 

3) 

(10/28/82) 

6 

6 

Normal 

01B 

1 

7 

(2) 

2 

8 

3 

9 

4 

10 

5 

11 

b 

12 

02B 

t 

13 

(3) 

2 

14 

3 

15 

4 

16 

5 

17 

6 

18 

02C 

1 

19 

(4) 

2 

20 

3 

21 

4 

22 

5 

23 

6 

24 

03  C 

1 

25 

(5) 

2 

26 

3 

27 

4 

28 

5 

29 

6 

30 

Note;  1)  Data  from  borehole  at  site  1  (60  a  hole)  to  SP  ch  02  (in  CTV). 

2)  Wind  speed  to  SP  ch  03  (in  CTV). 

3)  Wind  direction  to  SP  ch  05  (in  CTV). 


TABLE  IV. 1 

Status  of  NORSAR  SP  Instruments  recorded  on  data  tape. 
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St.  no. 

NS  (m) 

EW  (i 

2 

0 

0 

1 

271 

-62 

6 

-348 

-219 

9 

-407 

-660 

11 

1232 

-502 

TABLE  IV. 3 

Temporary  5-site  3-component  SP  seismometer 
experimental  array. 
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Fig.  IV. 1 


Task 


Date 


S A/Are a 


06C 

NORESS 

Activities  related  to  expansion 
of  the  NORESS  array 

( Apr. ) 

(Weeks  18,19  Ma 

06C 

NORESS 

Installation  of  3-component  SP- 
seism.  SS-1  on  five  new  stations 

(  "  20,21  Ma 

06C 

NORESS 

Activities  rel.  to  expansion  of 

NORESS  array  continued  involving 
attempt  to  find  a  possible  location 
for  the  RT  and  calculating  different 
fiber  optic  cable  lengths 

(June ) 

04C 

area 

Replaced  part  of  com.  cable 

(1  Jul) 

01 A 

Battery  bank  electrolyte  replaced 

(27  Jul) 

01 B 

_M  _ 

(  "  ) 

01A 

Noise  data  recording  (ch.  1  and 

60  m  hole) 

(start  28  Jul) 

02  B 

Battery  band  electrolyte  replaced 

(29  Jul) 

04  C 

Cable  splicing  SP00.  Protection  card 
replaced . 

(2  Aug) 

Q6C 

Battery  bank  electrolyte  replaced 

(3  Aug) 

01 B 

Cable  splicing  SP03  area 

(4  Aug) 

04C 

Cable  mending  near  SP01 

(8,9  Aug) 

01A 

RA-5  replacement 

(11  Aug) 

02B 

(12  Aug) 

04C 

RA-5  replacement  (03).  Battery 
bank  electrolyte  replaced 

(15  Aug) 

03C 

Battery  bank  electrolyte  replaced 

(16  Aug) 

04C 

Remaining  RA-5  ampl.  replaced 

(17  Aug) 

02C 

Battery  bank  electrolyte  replaced 

(18  Aug) 

04C 

Remote  Centering  Device  (RCD)  Vert. 

LP  seism,  replaced 

(22  Aug) 

06C 

Jobs  of  different  character  carried  ou 

t(25,2b  Aug) 

06C 

Different  tasks  related  to  profile 
blasting 

(29  Aug) 

02C 

RA-5  replacement  except  01,02 

(30  Aug) 

06C 

Profile  blasting  in  the  NORESS  array 
(5  men)  in  connection  with  a  seismic 
reflection  survey 

(30  Aug) 

NMC 

Property  control  by  J .  Guy  Turner 

Det  16  AFCMC,  Wiesbaden 

(12  Sept) 

03C 

SA  visit  in  connec.  with  power  failure 

(14  Sept) 

01B 

RA-5  ampl.  replaced  all  SP  points 

(20  Sept) 

03C 

-”-  (-SP04) 

(21  Sept) 

03C 

Faulty  cal.  amp.  cir.  repaired  (SPOl) 

(23  Sept) 

06C 

Reported  line  fault  checked 

(29  Sept) 

03C 

SA  visited  in  connection  with  loss  of 

(29  Sept) 

data;  found  power  line  broken 

TABLE  IV. 6 

Activities  in  the  field  and  at  NORSAR  Maintenance  Center. 
(1  April  -  30  September  1983) 


V.  DOCUMENTATION  DEVELOPED 

Doornbos,  D.J.:  )n  the  determination  of  radiated  energy  and  related  sour  e 
parameters,  submitted  for  publication. 

Husebye,  E.S.  and  E.  Thoresen:  Personal  seismometry  now,  submitted  for 
publication. 

Tronrud,  L.B.,  1983:  Semiannual  Technical  Summary,  NORSAR  Sci.  Rep. 

No.  2-82/83,  NTNK/NORSAR,  Kjeller,  Norway. 

Tryggvason,  K.  and  E.S.  Husebye:  Seismic  image  of  the  hypothesized  Icelandic 
hot  spot,  Tectonophysics ,  in  press. 


L.B.  Tronrud 
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VI.  SUMMARY  OF  TECHNICAL  REPORTS/ PAPERS  PREPARED 

VI. I  Example  of  a  seismic  absorption  band  at  high  frequencies 

There  is  growing  evidence  that,  in  most  regions  of  the  earth,  seismic 
absorption  in  terms  of  Q-^(u)  forms  a  band  centered  at  relatively  low 
frequencies,  i.e.,  the  high-frequency  cut-off  of  this  absorption  band 
is  in  the  range  of  short-period  body  waves  (0. 2-1.0  Hz)  (Doornbos,  1983). 
One  should  nevertheless  consider  the  possibility  that  the  absorption  band 
is  at  higher  frequencies  in  those  regions  which  are  known  to  be  highly 
absorptive,  or  'low-Q'  (e.g.,  Anderson  and  Given,  1982).  Examples  of 
such  regions  are  thought  to  be  the  upper  mantle  low-velocity  zone,  the 
region  near  the  base  of  the  mantle  (the  D"  layer)  and  the  region  near 
the  top  of  the  inner  core.  Observational  verification  is  usually  diffi¬ 
cult,  but  the  relative  amplitudes  and  waveforms  of  some  short-period 
core  phases  provide  an  interesting  case.  Short-period  NoRSAR  records 
show  a  characteristic  change  of  the  PKIKP  waveform  passing  through  the 
inner  core,  as  compared  to  PKPg£  bottoming  above.  One  can  attribute  this 
change  to  the  effect  of  absorption  in  the  inner  core,  provided  adequate 
care  is  taken  to  avoid  frequency-dependent  elastic  effects,  to  avoid 
source  effects  and  to  eliminate  receiver  structure  effects.  We  have 


used  full  NORSAR  array  beams  to  eliminate  near-receiver  effects,  and 
this  limits  our  analysis  to  events  between  1971-1976,  when  the  array 
still  had  its  full  size.  Fig.  VI. 1.1  shows  PKIKP  and  PKPg^  from  three 
such  events  in  the  distnce  range  197-151°.  The  two  phases  are  charac¬ 
teristically  different,  but  they  should  have  a  similar  waveform  apart 
from  absorption  effects.  By  verifying  that  the  Hilbert  transform  of  PKPgc 
is  similar  to  the  PKP^B  waveform,  it  can  be  concluded  that  the  anomalous 
waveshape  is  in  PKIKP.  To  investigate  this  effect  we  convolved  PKPgc 
with  an  absorption  band  operator  characterizing  transmission  through 
tlie  inner  core;  we  also  corrected  the  amplitude  ratio  due  to  purely 
elastic  effects.  The  required  inner  core  absorption  band  turns  out  to 
be  on  the  high-frequency  side  of  the  data.  Absorption  bands  encompassing 
the  data  (a  'constant  Q'  model)  or  on  the  low-frequency  side  of  it  (the 
more  usual  absorption  band)  are  unsuccessful.  Fig.  VI.1.2.a-c  illustrate 
the  effect  of  typical  examples  of  such  absorption  bands.  The  conclusion 
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is  that  for  absorption  in  the  inner  core,  the  low-frequency  cut-off  of 
the  band  is  above  2  Hz,  and  minimum  Qa  in  the  center  of  the 

band  is  likely  to  be  below  100.  The  intrinsic  dissipation  can  be  in  shear 
and/or  in  bulk.  Recently  proposed  bulk  dissipation  mechanisms  where  the 
governing  equations  are  diffusive  (Loper  and  Fearn,  1983)  have  consequences 
for  the  high-frequency  side  of  the  absorption  band.  However,  with  the  rela¬ 
tive  position  of  the  inner  core  absorption  band  as  inferred  here,  the  dif¬ 
ferences  are  unlikely  to  be  observable. 


D.J.  Doornbos 
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BC 


BC 


Fig.  VI. 1.1  NORSAR  array  beams  with  PKIKP(DF)  and  PKP(BC)  from  three 
events  in  the  distance  range  147-151°.  Relative  times 
adjusted  to  the  first  peak  or  trough.  Record  length  is 
5  seconds.  The  relative  amplitude  scaling  is  indicated 
to  the  right. 


w/I* 


OF  'BC' 


DF  ’BC’ 


a) 


Fig.  VI. 1.2  Observed  PKIKP(DF)  compared  with  synthetic  derived  from 

PKP(BC)  by  correcting  for  relative  path  effects,  including 
absorption  in  the  inner  core.  Other  details  as  in  Fig.  VI. 1.1. 
Inner  core  absorption  bands  are  specified  by.  (a)  t$=0.3  s, 
(2nr  )_I»0.0625  Hz,  ( 2 «r t ,  ) " 1  - 16  Hz;  (b)  t*=1.0  s, 
(2ht\)-1«0.002  Hz,  ( 2ntj )-1 "0 . 5  Hz;  (c)  t*-l.A  s, 

( 2 t!  c ^  ~ 1  =,i*  Hz’  Hz. 
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VI. 2  Spectral  bandwidth,  pulse  width  and  moments  in  source  analysis 
An  observational  parameter  often  employed  in  source  analysis  is  the 
spectral  bandwidth  (usually  inferred  from  a  corner  frequency  in  the 
spectrum);  for  example  radiated  energy  is  proportional  to  a  weighted 
average  of  the  spectral  bandwidth  squared.  On  the  other  hand,  a  measure 
of  source  size  is  given  by  the  time  domain  pulse  width  (more  precisely, 
the  second  central  moment  of  the  pulse).  Although  these  two  parameters 
are  simply  related  for  any  point  on  the  focal  sphere,  the  effect  of 
averaging  can  be  quite  different. 

It  is  possible  to  discuss  this  effect  in  the  context  of  a  moment 
tensor  representation  of  the  source.  If  the  temporal  derivative  of 
stress  glut  is  approximated  (Doornbos,  1982): 

mjk(£,T)  =  Mjk  fU,T)  (1) 

then  the  scalar  function  f(£,t)  can  be  expanded  in  moments.  If  we  choose 
as  the  reference  point  the  'center  of  gravity'  or  'centroid'  (£0>To^ 
of  f(£,x),  then  the  expansion  is  in  terms  of  the  central  moments.  We 
cancel  the  observationally  troublesome  phase  effect  by  relating  the 
moments  to  the  amplitude  density  spectrum  of  a  (normalized)  pulse 
f(i): 

1 

|f(w)  I2  -  1-W2  F(2)  +  —  u>“  (F(4)  +  3  F^2))  +  ...  (2) 

and  directivity  is  included  in  the  model  oy  replacing 

f(t)  by  f(T,£),  F(2)  by  £TF(2)C.  ,  etc. 

Here  5  is  a  generalized  slowness  vector  as  defined  in  Doornbos  (1982). 
For  practical  purposes  it  is  necessary  to  reduce  the  number  of  para¬ 
meters,  and  this  is  done  by  means  of  a  suitable  extrapolation  of  the 
spectrum  based  on  the  second  central  moments.  The  Gaussian  and  oo-square 
models  are  representative  examples. 
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The  corner  frequency  has  been  conventionally  determined  as  an  average 
over  the  focal  sphere.  Similarly  we  can  form  the  averaged  spectral 
bandwidth  lor  (’  and  S  waves: 


1 

Bc  =  —  ;  B(;f)  dft 

4  71  ii 


(3) 


where  c  is  the  P  or  S  velocity.  The  result  will  depend  on  the  direciivit 
of  f(£,t),  in  particular  on  the  effect  of  rupture  velocity.  In  contrast, 
the  averaged  pulse  width  squared  is  related  to 
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5c  F(2)  5c  dn  “ 
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+  Fzz> 


+  F 


tt 


(4) 


where  are  the  diagonal  components  of  F(2)>  This  result  does  not  ex¬ 
plicitly  involve  rupture  velocity.  As  a  measure  of  pulse  width  we  will 
use  L>c  =  n  F|,  with  x  a  (model  dependent)  constant.  Usually  an  inverse 

/N 

square  root  relation  between  B  and  F(2)  exists  for  any  point  on  the  focal 
sphere,  but  for  sources  with  strong  directivity  effects  (e.g.,  Haskell 
type  of  models),  the  result  of  averaging  procedures  (3)  and  (4)  can  be 
quite  different.  This  is  illustrated  by  comparing  Figs.  VI. 2.1  and  VI. 2. 2, 
which  give  results  for  moment  tensor  approximations  to  a  Haskell  bidirec¬ 
tional  model  with  aspect  ratio  0.4  and  to  a  circular  model,  respectively. 
The  'corner  frequency  shift'  in  terms  of  Ba/ Bg  is  also  illustrated 
in  these  figures  (c.f.  Hanks,  1981).  Fc  can  be  related  to  source  finite¬ 
ness,  and  Bc  to  the  dominant  frequency  range  for  energy  radiation. 


D.J.  Doornbos 
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DIRECTIVITY  COEFF.d 


Fig.  VI. 2.1  Results  for  averaged  spectral  bandwidth  and  pulse  width  Bc , 

Dc ,  c  the  P  or  S  velocity.  A  Gaussian  spectral  model  has  been 
used  in  the  moment  tensor  approximation  to  the  Haskell  bi¬ 
directional  model  with  aspect  ratio  0.4.  The  directivity 
coefficient  d  •  F^/(B  F£t).  For  any  point  on  the  focal 
sphere  2  Bc  Dc  *  1 ,  but  the  averaged  values  give  in  general 
2  Bc  Dc  >  1.  The  dotted  curve  denoted  <Ba/Bg>  gives  the 
ratio  of  areas  of  the  focal  sphere  where  Ba>Bg  and 
Ba<Bg,  respectively. 
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VI. 3  The  new  regional  array:  1983  field  experiments 

The  planning  work  for  the  new  regional  array  to  he  installed  in  1984  has 
continued  during  the  reporting  period.  There  have  been  extensive  dis¬ 
cussions  between  NORSAR  personnel  and  representatives  of  the  Sandia 
Laboratories,  Albuquerque,  New  Mexico,  which  is  the  U.S.  organization 
responsible  for  supplying  and  installing  all  hardware  (seismometers, 
amplifiers,  power  cables,  fiber  optic  cables  for  signal  transmission, 
electronic  equipment  at  the  central  site)  for  the  new  array.  The  new 
array  will  be  depLoyed  around  NORSAR  station  06C02  and  according  to 
status  of  the  project  as  of  November  1983  the  new  array  should  be 
operational  by  September  1984.  This  contribution  gives  a  description 
of  experiments  with  temporary  field  installations  during  the  summer 
of  1983  and  some  preliminary  results  from  analysis  of  the  data  col¬ 
lected.  These  experiments  were  undertaken  in  order  to  have  a  'last 
minute'  check  on  current  design  ideas  for  the  1984  array. 

Bore  ho  le_exper  iment s 

Results  from  a  study  of  high-frequency  noise  recorded  at  the  bottom  of 
a  60  ra  deep  borehole  and  comparison  with  simultaneously  recorded  noise 
at  the  surface  have  previously  been  reported  by  Bungum  (1983).  During 
1983,  PDR-2  recording  equipment  was  operated  in  trigger  mode  to  collect 
event  data  for  the  same  experimental  configuration.  An  example  showing  the 
P-phase  from  a  local  event  at  a  distance  of  about  3°  is  given  in  Fig. 

VI. 3.1.  As  is  seen  from  the  scaling  factors  the  maximum  amplitude  of 
the  signal  Is  drastically  reduced  for  the  borehole  record.  Spectral 
differences  are  particularly  pronounced  cbove  8-10  Hz.  It  is  suggested 
that  the  signal  loss  in  the  borehole  is  due  to  destructive  interference 
of  the  direct  P-wave  with  the  surface  reflected  one.  Crude  calculations 
on  wavelengths  and  time  delays  involved  tend  to  support  this  assumption. 

The  problem  at  hand  is  that  of  finding  the  optimum  depth  of  deployment 
ot  the  3-axis  package  to  go  Into  a  borehole  at  the  center  of  the  new 
array.  Data  collected  as  described  above  in  addition  to  borehole  data 
collected  in  the  U.S.  will  be  evaluated  to  settle  this  question.  An 
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uncertainL  uotor  Her.-  is  Che  coupling  between  the  seisaoiueters  and 
tlie  bednck  mil  its  influence  on  recorded  amplitudes  at  different 
f  requeue i o  s . 

3-axls  expe  r  iment 

During  the  period  June  LO-July  5,  1983,  data  were  recorded  on  a  5- 
statlon  array  of  3-axis  stations.  The  geometry  of  that  array  is  sho.  l 
in  Fig.  TV.  1  .  During  this  period  a  number  of  local  events  were  recorded. 

The  new  array  to  be  installed  in  1984  will  comprise  4  3-axis  elements, 
whereof  one  will  be  located  in  a  borehole  at  the  center  of  the  array. 

The  location  of  the  remaining  3  3-axis  elements  is  still  under  considera¬ 
tion,  and  any  one  of  the  24  remaining  sites  is  and  remains  a  possible  can¬ 
didate  for  deployment  of  3-axis  systems.  The  data  collected  this  summer 
will  now  be  subjected  to  analysis  with  the  purpose  of  deciding  where 
to  deploy  the  3  3-axis  sets  In  the  new  array  next  summer. 

2_l^channel  array  experiment 

Theoretical  work,  related  to  the  design  of  the  25-eleraent  vertical  arrav 
has  been  discussed  by  Mykkeltveit  11983).  Discussions  were  held  during 
the  spring  and  early  summer  of  L983  between  NORSAR  personnel  and  repre¬ 
sentatives  >f  Lawrence  Livermore  National  Laboratory  on  the  subject  of 
the  geometry  of  the  new  array.  It  was  agreed  to  test  the  proposed 
geometry  by  a  temporary  Installation  during  the  summer  of  1984.  The 
temporary  array  became  operational  on  July  23,  1983,  and  will  he  operated 
throughout  the  winter  of  1984.  It  comprises  21  channels  and  its  geometry 
is  identical  to  that  of  the  proposed  25-element  array  with  the  exception 
of  4  channels  in  the  outer  ring  which  have  not  been  installed.  The  geo¬ 
metry  of  the  temporary  array  is  shown  in  Fig.  IV. 2.  Analysis  of  the 
data  recorded  so  far  indicates  that  the  geometry  agreed  upon  is  a  very 
useful  one,  and  the  final  decision  has  been  made  to  Implement  it  in  1984. 

The  correlation  curves  for  signals  and  noise  on  which  the  proposed  geo¬ 
metry  was  base.!  Have  been  confirmed  by  the  new  data  as  shown  in  Fig. 
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VI. 3. 2,  providing  a  much  denser  sampling  and  broader  range  of  intersensor 
spacings  than  available  at  earlier  times. 

The  performance  of  the  f-k  analysis  for  detected  phase  arrivals  will 
be  of  crucial  importance  to  the  on-line  location  capability  of  the  new 
array.  We  are  now  in  the  process  of  collecting  and  analyzing  data  from 
events  with  known  location  (or  location  known  to  within  a  few  km). 

Results  of  the  f-k  analysis  (azimuth  and  phase  velocity)  for  regional 
events  analyzed  so  far  are  very  promising.  An  example  is  given  in 
Fig.  VI. 3. 3  showing  results  from  f-k  analysis  of  the  first  P  arrival 
(Pg)  from  a  local  event.  F-k  analysis  is  performed  at  b  and  8  Hz,  and 
the  true  azimuth  to  the  event  is  200°,  which  is  also  the  azimuth  re¬ 
sulting  fro n  tile  f-k  analysis  for  both  frequencies.  For  more  results 
here,  see  Section  VI. 4. 


S.  Mykkeltveit 
H.  Bungun 
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Fig*  VI. 1.2  Noise  correlation  curves  for  two  frequency  bands  for  the  21-element  array. 

The  data  interval  is  102*4  sec  long.  Mean  values  and  standard  deviations 
within  separation  intervals  of  ISO  m  are  y,iven  by  special  symbols. 


of  200°  and  distance  118  km  subjected  to  f-k  analysis.  The  data  uindou  (between 
vertical  bars)  for  f-k  analysis  contains  the  onset  of  the  f'g  phase.  Data  have 
been  filtered  3-9  Hz  before  analysis. 


VI. 4  Further  RONAPP  developments 

The  RONAPP  -  Regional  On-Line  Array  Processing  Package  -  has  now  been 
developed  further,  and  with  two  main  changes: 

1.  To  analyze  the  data  from  the  new  21-channel  test  array  in  the  NORESS 
siting  area  (see  Section  VI. 5) 

2.  To  include  beam  forming  in  the  detector. 

The  first  of  these  tasks  was  a  trivial  software  job,  but  the  second  was  not. 
To  summarize  what  has  been  reported  on  previously  (Mykkeltveit  et  al,  1982; 
Mykkeltvelt  and  Bungum,  1983),  the  basic  logic  of  RONAPP  is  as  follows: 

1)  Initialize  program,  parameters,  etc. 

2)  Event  j)e£e£t£0£  P.r.2c.£dur£:  Enter  a  Detection  Processor  and  stay  there 
until  a  detection  is  found 

3)  Event  J)e£e£tjLo£  Anal^sl_s:  When  a  detection  is  found,  analyze  the 
wave  train  causing  the  detection,  essentially  with  the  purpose  of 
finding  arrival  azimuth  and  phase  velocity 

4)  Event  !£ca£t£n_p£0£ed_ure :  If  the  wave  train  can  be  identified  as 
an  S  (or  Lg)  phase,  check  with  previous  detections  for  a  matching  P 
phase,  and  perform  an  event  location  if  possible 

5)  Return  to  2)  and  continue. 

Up  to  now,  most  of  the  RONAPP  tests  have  been  performed  with  data  from 
a  6-channel  test  array,  and  an  effort  at  that  stage  to  include  beam  forming 
in  the  detector  was  not  particularly  successful  because  of  beam  space 
instabilities  and  associated  problems  with  detection  reductions.  Because 
of  this  effect  of  a  poor  array  configuration,  most  of  the  initial  RONAPP 
analysis  was  based  on  detections  from  on1 y  one  vertical  beam. 

The  more  recent  changes  in  RONAPP  can  be  described  as  follows: 

1)  The  detection  processor  (DP)  is  initialized  with  any  number  of  beams, 
each  one  specified  in  terms  of  azimuth,  inverse  velocity,  filter,  and 
individual  channel  weights.  Time  delays  are  then  computed  once  and 
for  all,  and  it  will  be  easy  at  a  later  stage  to  include  possible 
time  delay  corrections.  The  filter  is  normally  the  same  for  all  beams 
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because  of  the  beam  space  instabilities  that  otherwise  would  occur. 

The  use  of  several  filters  would  require  equally  many  independent 
DP  partitions,  and  this  would  only  be  a  question  of  computer  time. 

2)  The  STA  computations  are  now  based  on  averaging  absolute  instead  ot 
squared  amplitudes,  in  order  to  avoid  precision  problems  in  the  computer 
and  to  save  computer  time.  The  corresponding  loss  In  SNR  is  negligible. 

3)  For  each  beam,  there  is  one  STA/LTA  threshold  for  declaring  a  detection 
(provided  that  P  out  of  Q  successive  samples  exceeds  the  threshold) 

and  another  (and  lower)  threshold  for  closing  the  detection. 

4)  The  DP  is  defined  to  be  In  detection  state  if  a  detection  is  declared 
for  at  least  one  beam.  To  leave  the  detection  state  (i.e.,  to  a '  low 

a  new  detection  to  occur),  two  criteria  must  be  fulfilled: 

All  but  a  specified  number  of  beams  must  be  out  of  their  (indi¬ 
vidual)  detection  state.  (That  number  has  so  far  mostly  bee  <  set 
to  zero.) 

A  certain  time  must  have  elapsed  since  the  last  beam  closed  its 
detection. 

■j)  To  process  a  detection  (i.e.,  to  enter  Event  Detection  Analysis), 

a  certain  time  must  have  elapsed  since  the  previous  processed  detection. 

6)  Before  a  detection  is  processed,  the  beam  with  maximum  STA  is  found 
and  used  as  a  basis  for  determination  of  refined  arrival  time,  dominant 
frequency,  f-k  analysis  prefilter,  and  f-k  analysis  time  window. 

7)  Following  each  f-k  analysis,  the  RONAPP  procedure  is  unchanged,  i.e., 
a  phase  association  and  event  location  procedure  is  entered  if  the 
last  detection  has  been  Identified  as  an  S-type  phase.  So  far,  our  data 
base  is  too  small  for  development  of  possible  regional  corrections  In 
locations  (systematic  azimuth  and  phase  velocity  deviations),  but  with 
the  present  program  structure  the  inclusion  of  such  corrections  shou  d 
be  quite  straightforward  software-wise. 

This  version  of  the  RONAPP  package  has  been  tested  on  a  number  of  selected 

events,  as  well  as  on  real-time  data.  The  performance  with  respect  to  an 

event  (explosion)  in  western  Norway  is  demonstrated  in  the  following,  with 

the  seismic  data  shown  in  Fig.  VI. 4.1,  and  detections  indicated  by  .  rrow;  . 
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We  have  previously  had  some  problems  with  too  many  o»la  detections,  a 
situation  which  now  has  been  significantly  improved  with  the  availability 
of  data  from  the  new  21-channel  array  and  with  the  new  detection  reduction 
procedures  discussed  above.  The  point  with  the  coda  detections  is  of  course 
that  they  should  be  reduced  in  number  without  losing  the  detections  that 
we  are  interested  in. 

The  results  for  each  detection  are  shown  in  the  detection  report  on  top 
of  Table  VI. 4.1,  and  at  the  bottom  of  that  table  the  results  from  phase 
association  and  location  are  given.  It  is  seen  there  tha,.  the  first  location 
is  based  on  the  Sn  phase  combined  with  the  Pn  phase,  a  location  which  is 
recomputed  upon  the  arrival  of  the  stronger  Lg  phase.  This  is  exactly 
the  way  we  want  the  processing  package  to  perform. 

In  Table  VI. 4. 2  the  results  tied  to  one  particular  detection,  namely, 
the  Lg  phase,  are  demonstrated.  There  are  11  beams  (1  is  vertical,  2-7  are 
P  beams,  8-13  are  Lg  beams)  which  all  have  detected,  but  it  is  interesting 
to  note  that  the  best  beam  (no.  12)  has  an  SNR  of  140  while  many  of  the 
others  have  around  30.  This  gives  a  ratio  of  about  4.7  which  corresponds 
almost  exactly  to  the  beamforming  gain  that  should  be  expected  under  ideal 
signal  and  noise  conditions.  It  is  seen  further  down  in  Table  VI. 4. 2 
that  beam  12  Is  located  (in  slowness  space)  very  close  to  where  the  f-k 
analysis  finds  maximum  power. 


H.  Bungum 
S.  Mykkeltveit 
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Table  VI. 4. 2 


Detector  output  for  the  last  (Lg)  detection  in  Figs.  VI. 4. 1-2. 
Beam  12.  with  largest  SNR,  has  an  azimuth  of  240°  and  a  phase 
velocity  of  4.5  kra/s. 
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VI. 5.  Weighted  beamforming  In  a  real  time  environment 
As  is  well  known,  the  extent  of  signal  and  noise  correlation  between 
sensors  in  an  array  might  significantly  affect  its  performance  in  terms 
of  suppressing  the  ambient  noise  while  retaining  signal  integrity.  The 
importance  of  signal  and  noise  correlations  as  a  function  of  sensor 
separations  and  frequency  band  has  recently  been  demonstrated  by 
Mykkeltveit  et  al  (1983)  in  a  scheme  for  optimizing  array  configurations. 
A  rather  obvious  result  here  is  that  once  the  array  becomes  operational 
(configuration  fixed)  its  performance  is  much  dependent  on  dominant  sig- 
gnal  frequency.  For  example,  a  miniarray  designed  for  optimum  detection 
capabilities  for  signals  from  local  and  regional  events,  say  in  the  range 
5-15  Hz,  would  be  far  less  efficient  for  teleseismic  signals  in  the  fre¬ 
quency  range  1-3  Hz.  From  a  general  seisaological  point  of  view,  the 
preference  is  for  array  operation  which  is  not  strongly  peaked  as  a 
function  of  frequency.  In  practice  this  requires  that  a  flexible  weight¬ 
ing  scheme  is  introduced  as  part  of  the  array  on-line  operation,  where 
we  try  to  capitalize  on  the  information  contained  in  the  noise  covariance 
matrix.  The  weighting  technique  used,  to  be  briefly  described  in  the 
following,  is  very  similar  to  the  optimum  processing  schemes  developed 
by  Capon,  Lacoss  and  others  (e.g.,  see  Lacoss,  1974). 

S^£na_l^nolse_mod£lJLing_and_Oj3t:iraura_we^£fa  estimation 

Many  weighting/filtering  schemes  have  been  developed  for  multichannel 
noise  suppression.  The  best  known  class  here  is  the  Wiener  filters 
which  utilize  the  information  contained  in  the  autocovariance  matrix, 
while  in  our  scheme  the  subset  hereof,  the  covariance  matrix,  is  in 
focus . 

Model  I  -  Signals  identical: 

y^  »  s  +  n|  ;  s  *  signal,  nj  »  noise  i-tl  sensor 

FX^nj)  -  T*j 


(T^j)  *  noise  covariance  matrix 
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Standard  (unweighted)  beams: 


§  =  4  =  ks  +  l  n^ 

i  =  l 

I  =  k2s2  +  i'll 


4l,k  = 
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Weighted  beam: 


§  =  wivi  =  s  V.  w^  +  E  w^nj: 
1  =  1 


Gain 


U'w)2 


w '  iw 


wik 


«1 


wk 


By  derivation  of  the  gain  function  we  find: 
max(Gain2)  for  w  ~ 


In  practice  we  want  to  use  a  normalized  weight  function,  and  this  intro¬ 
duces  the  wnorm-weights : 


«'w 


wuse  =  w,/wnorra 


Implicit  In  this  estimation  scheme  is  that  Iwj  =  1,  Important,  in 
order  to  avoid  negative  weights,  the  w^  elements  in  the  above  gain 
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function  are  replaced  by  «t  -  x^.  For  k  <  20  the  minimizing  of  the 
gain  function  (is  negative  counterpart)  car.  easily  be  performed  directly 
without  derivatives,  etc. 

Model  II  ~  Signals  not  perfectly  correlated 

In  the  case  the  gain  would  be: 
w'ftw 

Gain^  “  - 

w'Zw 

with  11  being  the  signal  covariance  function. 

Again,  we  can  show  that: 

max  Gain  for  (fl-AE)  w  “  0 

where  A  is  the  largest  eigenvalue  and  w  the  corresponding  eigenvector. 

We  may  also  extend  the  model  to  Include  a  priori  knowledge  of  signal 
amplitudes  (e.g.,  see  Christof fersson  and  Husebye,  1974).  Note  that 
significant  gain  in  beamforming  is  obtainable  when  sensor  no.se 
correlation  is  significant.  Restricting  weights  to  be  positive  only 
essentially  amounts  to  0/1  weighting. 

Practical  considerations 

As  in  all  other  noise  suppression  schemes,  the  critical  factor  is  noise 
stationarity,  in  this  particular  case  the  time/space  stability  of  the  noise 
covariance  estimate  in  use.  Parameters  of  importance  here  are: 

-  IW  “  length  of  data  window  for  covariance  estimate  (2.5-10  sec); 

in  average  most  stable  results  for  IW  *  10  sec  or  400  samples 

-  INCR  “  covariance  matrix  updating  rate;  typical  values  1.25- 

2.50  sec;  in  average  most  stable  results  for  INCR  “  2.5  sec 
or  100  samples 

-  ALFA  “  covariance  matrix  smoothing  parameter;  similar  functioning 

of  the  STA/LTA  sliding  windows  as  used  for  signal  detectors; 
here  defined  as  ALFA  -  1.  -  EXP(AL0G(1.-PR0)/FIW)  with 
PRO  -  0.5  and  FIW  -  IW  +  1 
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Now,  parameters  of  importance  in  judging  the  performance  of  the  above 

weighting  scheme  are  as  follows: 

TN  -  theoretical  gain  permitting  negative  weights  tied  to  the 
covariance  matrix  only;  no  real  data  used 
ON  -  optimum  gain  (neg.  weights  permitted)  as  measured  on  the 

noise 

TP  =  theoretical  gain  for  positive  weights  only 

-  pw  =  optimum  gain  for  positive  weights  only 

-  ZW  “  only  positive  weights  larger  than  0.05  Included 

-  SB  “  standard  beam  gain  in  slgnal-to-noise  ratio  (SNR). 

The  above  gain  measures  are  given  in  dB  relative  SB.  It  should  be  remembered 
that  it  is  important  to  mask  properly  'dead*  traces,  as  these  otherwise 
would  be  given  very  large  weights  in  view  of  their  vanishing  variances. 

On  the  other  hand,  spikes  would  be  efficiently  removed  as  such  traces  have 
large  variances. 


Real-time  simulation_of_  weijghtinfc  scheme 

To  test  the  performance  of  the  above  weighting  scheme,  a  limited  number 
of  noise  samples  for  the  new.  prototype  NORESS  array  (see  Fig.  VI. 5.1)  have 
been  used  in  data  analysis.  Due  to  disk  storage  restriction,  only  13 
channels  have  been  used  at  any  type:  1)  sensor  1,  B  &  C  rings  and  ii) 
sensor  1,  C  &  D  rings.  The  Inner  A  ring  sensors  have  been  consistently 
deleted  due  to  small  ranges,  as  their  sensor  separations  are  relatively 
small  vis-a-vis  the  upper  frequency  'cut-off'  at  about  4  Hi.  These  sensors 
may  well  be  used,  but  presently  strong  transmission  (electronic)  noise 
prevails  above  this  frequency,  i.e.,  the  provisional  60  dB  dynamic  range 
recording  system  does  not  permit  adequate  noise  sampling  over  a  wide 
frequency  range.  Typical  noise  data  used  in  analysis  is  displayed  in 
Fig.  VI. 5. 2;  filter  setting  1. 0-3.0  Hz  3rd  order  Butterworth.  The  strong 
correlation  between  sensors  in  the  Innermost  rings  is  rather  obvious. 

NoiS£  suppressiG n  capabilities 

The  optimum  gain  function  (negative  weights)  is  visually  displayed 
in  Fig.  VI. 5. 2,  while  more  detailed  results  are  presented  in  Tab' 


Vi. 5.1  for  a  typical  noise  sample.  Not  unexpectedly,  the  largest  noise 
suppression  is  obtained  for  low  frequencies,  that  is,  where  the  noise 
is  fairly  well  correlated.  This  is  also  shown  in  the  rightmost  column 
in  the  table,  where  the  number  of  times  a  sensor  has  contributed  con¬ 
structively  to  the  weighted  beamforming  is  indicated.  This  is  tied  to 
the  'positive'  weighting  scheme,  but  the  weights  contributing  by  less 
than  52  are  zeroed  while  the  remaining  are  set  equal  to  1.  The  total 
number  of  updating  intervals  was  44,  so  some  sensors  are  truly  abundant 
The  corresponding  relative  gain  is  listed  in  the  dW— column. 


Some  interesting  features  of  Table  VI. 5.1  are  as  follows: 
the  noise  suppression  is  azimuth-dependent 

the  weighting  gain  appears  to  be  velocity-dependent,  relatively 
best  performances  for  surface  wave  phase  velocities  (4.6  km  s  * 
and  teleseismic  P-velocities  (12  km  s_1) 

-  the  weighting  scheme  can  produce  gains  in  excess  of  <TT. 


A  similar  experiment  was  performed  by  replacing  the  Bring  sensors 
with  those  of  the  D-ring  (see  Table  VI. 5. 2).  The  outstanding  feature 
here  is  the  gain  from  weighted  beamforming  is  very  small,  occasionally 
even  negative.  The  reason  for  this  is  non-stationarity  in  the  noise  fiel. 
when  sensor  separations  become  relatively  large.  The  theoretical  gain  in 
this  case  is  not  too  different  from  that  of  the  previous  and  furthermore 
well  above  1  dB.  The  latter  data  are  not  shown  here. 


The  noise  suppression  is  from  Tables  VI. 5.1  6  VI. 5.2  o'  iously  azimuth- 
dependent,  and  this  feature  has  been  further  lnvestiga' ed.  Im  Table  VI. 5.) 
gala  as  a  function  of  small  azimuth  intervals  is  iispi  -yed.  Interestingly, 
for  those  azimuths  where  the  gain  (SB)  is  the  sma.lest,  the  most  contri¬ 
buting  sensors  (underlined  in  the  rightmost  coins)  for  triangles  .square 
whose  longest  axis  are  roughly  perpendicular  to  the  azi  uth  in  question 
(the  OK  notation  in  the  AC  coliam).  This  we  Interpret  i  terns  of  stronger 
correlation  in  the  noise  in  the  direction  of  propagation  (radial)  than 
in  the  transverse  direction,  which  Is  Indeed  expeted  f.r  scattering 


phenomena  (Chernov,  1960).  In  other  words,  the  correlation  In  the 
noise  appears  to  depend  on  both  sensor  separation  and  their  relative 
orientation.  Th.  latter  feature,  If  prominent,  would  clearly  have  a 
bearing  on  the  array  configuration  optimization  schemes.  For  example, 
the  approach  of  Mykkeltveit  et  al  (1983)  essentially  amounts  to  a  noise 
covariance  matrix  model  with  zeros  or  small  negative  values  for  all  off- 
diagonal  elements.  We  have  not  observed  this  during  our  analysis  so 
far  nor  obtained  a  gain  exceeding  /N  for  standard  beamforming  (SB). 

This  Is  a  bit  puzzling  as  noise  correlations  as  a  function  of  sensor 
separations  only  have  been  observed  to  exhibit  a  small  negative  minimum 
of  the  order  of  0.1  units.  This  discrepancy  may  reflect  a  certain  non- 
stationarity  in  the  noise  field  in  short  intervals  as  used  here,  although 
the  above-ment tuned  noise  covariance  matrix  model  is  likely  to  represent 
an  oversimplification  of  the  real  noise  field.  Thus,  if  the  off-diagonal 
elements  cannot  be  ensured  to  be  zero  or  negative,  the  consequence  of 
this  is  that  the  center  station  seldom  would  contribute  constructively 
towards  noise  suppression  by  beamforming.  This  is  quite  obvious  from 
the  results  presented  in  Tables  VI. 5.1-3  for  the  lower  frequency  bands. 

Dijscussi  on_and_con£l  us  ion 

A  preliminary  evaluation  of  the  provisonal  regional  array  installed 
within  the  NORSAR  array  has  been  *ade.  It  is  somewhat  incomplete  in 
the  sense  that  the  operational  settings  of  the  array  up  to  now  only 
permit  analysis  In  the  frequency  band  1  to  It  Hz  and  thus  excludes  the 
most  important  frequency  banc  fo.  small  local  <nd  regional  events, 
say  1.5  to  3.5  mb  units,  namely,  '  to  15  Hz.  Anyaay,  the  main  results 
obtained  are  as  follows: 

When  there  Is  a  significant  tie  in  the  noise  data,  weighted 
beamforming  would  produce  an  additional  noise  suppression  gain 
of  the  order  of  1  to  3  dB,  roughly  equivalent  to  10-30X  operational 
improvement. 

The  noise  field  must  be  rated  non-statlonary  as  significant  changes 
in  the  noise  covariance  take  place  within  5  to  10  sec  interval. 

In  this  respect  claimed  long  term  noise  stationarity  features 
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like  negative  correlation  values  in  certain  sensor  separation 
ranges  appear  to  be  very  modest,  at  least  in  comparison  to  rela¬ 
tively  strong  short  term  variations. 

A  certain  directabllity  in  the  noise  field  is  apparent  for 
certain  azimuth  directions.  In  practice  this  gives  reduced 
noise  suppression  capabilities  or  an  equivalent  higher  false 
alarm  rate  on  certain  beams. 

Array  configuration  optimization;  this  is  a  problem  if  reasonable 
performance  is  desired  over  a  relatively  wide  frequency  range. 
Optimum  processing  schemes  appear  to  be  unavoidable  here  or  the 
array  must  compromise  on  a  duality  in  configuration. 

Finally,  optimum  weighting  schemes  are  rather  time-consuming 
in  practice  and  can  hardly  be  'afforded'  unless  there  is  access 
to  an  array  processor  for  handling  the  covariance  matrix  'inversion' 
task.  We  are  pursuing  this  problem  now,  the  motivation  being  that 
this  would  be  very  cost  effective  per  unit  dB  gained  in  noise 
suppression  of  small  arrays. 


E.S.  Husebye 
S.K.  Ingate 

A.  Christof fersson,  Univ.  of 
Uppsala 
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